Lenalidomide Stabilizes Protein-Protein Complexes by Turning Labile Intermolecular H-Bonds into Robust Interactions.

Targeted protein degradation is a promising therapeutic strategy, spearheaded by the anti-myeloma drugs lenalidomide and pomalidomide. These drugs stabilize very efficiently the complex between the E3 ligase Cereblon (CRBN) and several non-native client proteins (neo-substrates), including the transcription factors Ikaros and Aiolos and the enzyme Caseine Kinase 1α (CK1α,), resulting in their degradation. Although the structures for these complexes have been determined, there are no evident interactions that can account for the high efficiency of formation of the ternary complex. We show that lenalidomide's stabilization of the CRBN-CK1α complex is largely due to hydrophobic shielding of intermolecular hydrogen bonds. We also find a quantitative relationship between hydrogen bond robustness and binding affinities of the ternary complexes. These results pave the way to further understand cooperativity effects in drug-induced protein-protein complexes and could help in the design of improved molecular glues and more efficient protein degraders.

A fluorogenic probe for granzyme B enables in-biopsy evaluation and screening of response to anticancer immunotherapies.

Immunotherapy promotes the attack of cancer cells by the immune system; however, it is difficult to detect early responses before changes in tumor size occur. Here, we report the rational design of a fluorogenic peptide able to detect picomolar concentrations of active granzyme B as a biomarker of immune-mediated anticancer action. Through a series of chemical iterations and molecular dynamics simulations, we synthesize a library of FRET peptides and identify probe H5 with an optimal fit into granzyme B. We demonstrate that probe H5 enables the real-time detection of T cell-mediated anticancer activity in mouse tumors and in tumors from lung cancer patients. Furthermore, we show image-based phenotypic screens, which reveal that the AKT kinase inhibitor AZD5363 shows immune-mediated anticancer activity. The reactivity of probe H5 may enable the monitoring of early responses to anticancer treatments using tissue biopsies.

Discovery and In Vivo Proof of Concept of a Highly Potent Dual Inhibitor of Soluble Epoxide Hydrolase and Acetylcholinesterase for the Treatment of Alzheimer's Disease.

With innumerable clinical failures of target-specific drug candidates for multifactorial diseases, such as Alzheimer's disease (AD), which remains inefficiently treated, the advent of multitarget drug discovery has brought a new breath of hope. Here, we disclose a class of 6-chlorotacrine (huprine)-TPPU hybrids as dual inhibitors of the enzymes soluble epoxide hydrolase (sEH) and acetylcholinesterase (AChE), a multitarget profile to provide cumulative effects against neuroinflammation and memory impairment. Computational studies confirmed the gorge-wide occupancy of both enzymes, from the main site to a secondary site, including a so far non-described AChE cryptic pocket. The lead compound displayed in vitro dual nanomolar potencies, adequate brain permeability, aqueous solubility, human microsomal stability, lack of neurotoxicity, and it rescued memory, synaptic plasticity, and neuroinflammation in an AD mouse model, after low dose chronic oral administration.

Dual Effect of Prussian Blue Nanoparticles on Aß40 Aggregation: ß-Sheet Fibril Reduction and Copper Dyshomeostasis Regulation.

Alzheimer's disease (AD), affecting almost 50 million individuals worldwide, is currently the first cause of dementia. Despite the tremendous research efforts in the last decade, only four supportive or palliative drugs, namely, acetylcholinesterase (AChE) inhibitors donepezil, galantamine, and rivastigmine and the glutamate NMDA receptor antagonist memantine, are currently available. New therapeutic strategies are becoming prominent, such as the direct inhibition of amyloid formation or the regulation of metal homeostasis. In the present report, the potential use of Prussian blue (PB), a drug that is in the World Health Organization Model List of Essential Medicines, in AD treatment is demonstrated. Both in vitro and in cellulo studies indeed suggest that PB nanoparticles (PBNPs) are capable of reducing the formation of typical amyloid-ß fibers (detected by thioflavin T fluorescence) and restoring the usual amyloid fibrillation pathway via chelation/sequestration of copper, which is found in high concentrations in senile plaques.

Combining Virtual Reality Visualization with Ensemble Molecular Dynamics to Study Complex Protein Conformational Changes.

Molecular dynamics (MD) simulations are increasingly used to elucidate relationships between protein structure, dynamics, and their biological function. Currently, it is extremely challenging to perform MD simulations of large-scale structural rearrangements in proteins that occur on millisecond timescales or beyond, as this requires very significant computational resources, or the use of cumbersome "collective variable" enhanced sampling protocols. Here, we describe a framework that combines ensemble MD simulations and virtual reality visualization (eMD-VR) to enable users to interactively generate realistic descriptions of large amplitude, millisecond timescale protein conformational changes in proteins. Detailed tests demonstrate that eMD-VR substantially decreases the computational cost of folding simulations of a WW domain, without the need to define collective variables a priori. We further show that eMD-VR generated pathways can be combined with Markov state models to describe the thermodynamics and kinetics of large-scale loop motions in the enzyme cyclophilin A. Our results suggest eMD-VR is a powerful tool for exploring protein energy landscapes in bioengineering efforts.

Investigating Cryptic Binding Sites by Molecular Dynamics Simulations.

This Account highlights recent advances and discusses major challenges in investigations of cryptic (hidden) binding sites by molecular simulations. Cryptic binding sites are not visible in protein targets crystallized without a ligand and only become visible crystallographically upon binding events. These sites have been shown to be druggable and might provide a rare opportunity to target difficult proteins. However, due to their hidden nature, they are difficult to find through experimental screening. Computational methods based on atomistic molecular simulations remain one of the best approaches to identify and characterize cryptic binding sites. However, not all methods are equally efficient. Some are more apt at quickly probing protein dynamics but do not provide thermodynamic or druggability information, while others that are able to provide such data are demanding in terms of time and resources. Here, we review the recent contributions of mixed-solvent simulations, metadynamics, Markov state models, and other enhanced sampling methods to the field of cryptic site identification and characterization. We discuss how these methods were able to provide precious information on the nature of the site opening mechanisms, to predict previously unknown sites which were used to design new ligands, and to compute the free energy landscapes and kinetics associated with the opening of the sites and the binding of the ligands. We highlight the potential and the importance of such predictions in drug discovery, especially for difficult ("undruggable") targets. We also discuss the major challenges in the field and their possible solutions.

Dynamic design: manipulation of millisecond timescale motions on the energy landscape of cyclophilin A.

Proteins need to interconvert between many conformations in order to function, many of which are formed transiently, and sparsely populated. Particularly when the lifetimes of these states approach the millisecond timescale, identifying the relevant structures and the mechanism by which they interconvert remains a tremendous challenge. Here we introduce a novel combination of accelerated MD (aMD) simulations and Markov state modelling (MSM) to explore these 'excited' conformational states. Applying this to the highly dynamic protein CypA, a protein involved in immune response and associated with HIV infection, we identify five principally populated conformational states and the atomistic mechanism by which they interconvert. A rational design strategy predicted that the mutant D66A should stabilise the minor conformations and substantially alter the dynamics, whereas the similar mutant H70A should leave the landscape broadly unchanged. These predictions are confirmed using CPMG and R1ρ solution state NMR measurements. By efficiently exploring functionally relevant, but sparsely populated conformations with millisecond lifetimes in silico, our aMD/MSM method has tremendous promise for the design of dynamic protein free energy landscapes for both protein engineering and drug discovery.

A computationally designed binding mode flip leads to a novel class of potent tri-vector cyclophilin inhibitors.

Cyclophilins (Cyps) are a major family of drug targets that are challenging to prosecute with small molecules because the shallow nature and high degree of conservation of the active site across human isoforms offers limited opportunities for potent and selective inhibition. Herein a computational approach based on molecular dynamics simulations and free energy calculations was combined with biophysical assays and X-ray crystallography to explore a flip in the binding mode of a reported urea-based Cyp inhibitor. This approach enabled access to a distal pocket that is poorly conserved among key Cyp isoforms, and led to the discovery of a new family of sub-micromolar cell-active inhibitors that offer unprecedented opportunities for the development of next-generation drug therapies based on Cyp inhibition. The computational approach is applicable to a broad range of organic functional groups and could prove widely enabling in molecular design.

First homology model of Plasmodium falciparum glucose-6-phosphate dehydrogenase: Discovery of selective substrate analog-based inhibitors as novel antimalarial agents.

In Plasmodium falciparum the bifunctional enzyme glucose-6-phosphate dehydrogenase‒6-phosphogluconolactonase (PfG6PD‒6PGL) is involved in the catalysis of the first reaction of the pentose phosphate pathway. Since this enzyme has a key role in parasite development, its unique structure represents a potential target for the discovery of antimalarial drugs. Here we describe the first 3D structural model of the G6PD domain of PfG6PD‒6PGL. Compared to the human enzyme (hG6PD), the 3D model has enabled the identification of a key difference in the substrate-binding site, which involves the replacement of Arg365 in hG6PD by Asp750 in PfG6PD. In a prospective validation of the model, this critical change has been exploited to rationally design a novel family of substrate analog-based inhibitors that can display the necessary selectivity towards PfG6PD. A series of glucose derivatives featuring an α-methoxy group at the anomeric position and different side chains at position 6 bearing distinct basic functionalities has been synthesized, and their PfG6PD and hG6PD inhibitory activities and their toxicity against parasite and mammalian cells have been assessed. Several compounds displayed micromolar affinity (Ki up to 23 µM), favorable selectivity (up to > 26-fold), and low cytotoxicity. Phenotypic assays with P. falciparum cultures revealed high micromolar IC50 values, likely as a result of poor internalization of the compounds in the parasite cell. Overall, these results endorse confidence to the 3D model of PfG6PD, paving the way for the use of target-based drug design approaches in antimalarial drug discovery studies around this promising target.

Correction: Unveiling a novel transient druggable pocket in BACE-1 through molecular simulations: Conformational analysis and binding mode of multisite inhibitors.

[This corrects the article DOI: 10.1371/journal.pone.0177683.].

Tetrahydrobenzothiophene carboxamides: Beyond the kinase domain and into the fatty acid realm.

A series of tetrahydrobenzothiophene carboxamides, inspired by structural features present in kinase and SCD1 inhibitors, are presented here. Prototype compound 8 (MMDD13) modulates fatty acid elongase and desaturase indexes, lipid accumulation, while preserving kinase inhibitory activity. This chemotype represents a stepping stone towards chemical probes to study the consequences of lipid metabolism modulation through non-redundant pathways.

Unveiling a novel transient druggable pocket in BACE-1 through molecular simulations: Conformational analysis and binding mode of multisite inhibitors.

The critical role of BACE-1 in the formation of neurotoxic ß-amyloid peptides in the brain makes it an attractive target for an efficacious treatment of Alzheimer's disease. However, the development of clinically useful BACE-1 inhibitors has proven to be extremely challenging. In this study we examine the binding mode of a novel potent inhibitor (compound 1, with IC50 80 nM) designed by synergistic combination of two fragments-huprine and rhein-that individually are endowed with very low activity against BACE-1. Examination of crystal structures reveals no appropriate binding site large enough to accommodate 1. Therefore we have examined the conformational flexibility of BACE-1 through extended molecular dynamics simulations, paying attention to the highly flexible region shaped by loops 8-14, 154-169 and 307-318. The analysis of the protein dynamics, together with studies of pocket druggability, has allowed us to detect the transient formation of a secondary binding site, which contains Arg307 as a key residue for the interaction with small molecules, at the edge of the catalytic cleft. The formation of this druggable "floppy" pocket would enable the binding of multisite inhibitors targeting both catalytic and secondary sites. Molecular dynamics simulations of BACE-1 bound to huprine-rhein hybrid compounds support the feasibility of this hypothesis. The results provide a basis to explain the high inhibitory potency of the two enantiomeric forms of 1, together with the large dependence on the length of the oligomethylenic linker. Furthermore, the multisite hypothesis has allowed us to rationalize the inhibitory potency of a series of tacrine-chromene hybrid compounds, specifically regarding the apparent lack of sensitivity of the inhibition constant to the chemical modifications introduced in the chromene unit. Overall, these findings pave the way for the exploration of novel functionalities in the design of optimized BACE-1 multisite inhibitors.

Mechanism of the Pseudoirreversible Binding of Amantadine to the M2 Proton Channel.

The M2 proton channel of influenza A virus is an integral membrane protein involved in the acidification of the viral interior, a step necessary for the release of the viral genetic material and replication of new virions. The aim of this study is to explore the mechanism of drug (un)binding to the M2 channel in order to gain insight into the structural and energetic features relevant for the development of novel inhibitors. To this end, we have investigated the binding of amantadine (Amt) to the wild type (wt) M2 channel and its V27A variant using multiple independent molecular dynamics simulations, exploratory conventional metadynamics, and multiple-walkers well-tempered metadynamics calculations. The results allow us to propose a sequential mechanism for the (un)binding of Amt to the wt M2 channel, which involves the adoption of a transiently populated intermediate (up state) leading to the thermodynamically favored down binding mode in the channel pore. Furthermore, they suggest that chloride anions play a relevant role in stabilizing the down binding mode of Amt to the wt channel, giving rise to a kinetic trapping that explains the experimentally observed pseudoirreversible inhibition of the wt channel by Amt. We propose that this trapping mechanism underlies the inhibitory activity of potent M2 channel blockers, as supported by the experimental confirmation of the irreversible binding of a pyrrolidine analogue from electrophysiological current assays. Finally, the results reveal that the thermodynamics and kinetics of Amt (un)binding is very sensitive to the V27A mutation, providing a quantitative rationale to the drastic decrease in inhibitory potency against the V27A variant. Overall, these findings pave the way to explore the inhibitory activity of Amt-related analogues in mutated M2 channel variants, providing guidelines for the design of novel inhibitors against resistant virus strains.

Blinded predictions of binding modes and energies of HSP90-α ligands for the 2015 D3R grand challenge.

In the framework of the 2015 D3R inaugural grand challenge, blind binding pose and affinity predictions were performed for a set of 180 ligands of the Heat Shock Protein HSP90-α protein, a relevant cancer target. Spectral clustering was used to rapidly identify alternative binding site conformations in publicly available crystallographic HSP90-α structures. Subsequently, multiple docking and scoring protocols employing the software Autodock Vina and rDock were applied to predict binding modes and rank order ligands. Alchemical free energy calculations were performed with the software FESetup and Sire/OpenMM to predict binding affinities for three congeneric series subsets. Some of the protocols used here were ranked among the top submissions according to most of the evaluation metrics. Docking performance was excellent, but the scoring results were disappointing. A critical assessment of the results is reported, as well as suggestions for future similar competitions.

Design, synthesis and biological evaluation of N-methyl-N-[(1,2,3-triazol-4-yl)alkyl]propargylamines as novel monoamine oxidase B inhibitors.

Different azides and alkynes have been coupled via Cu-catalyzed 1,3-dipolar Huisgen cycloaddition to afford a novel family of N1- and C5-substituted 1,2,3-triazole derivatives that feature the propargylamine group typical of irreversible MAO-B inhibitors at the C4-side chain of the triazole ring. All the synthesized compounds were evaluated against human MAO-A and MAO-B. Structure-activity relationships and molecular modeling were utilized to gain insight into the structural and chemical features that enhance the binding affinity and selectivity between the two enzyme isoforms. Several lead compounds, in terms of potency (submicromolar to low micromolar range), MAO-B selective recognition, and brain permeability, were identified. One of these leads (MAO-B IC50 of 3.54µM, selectivity MAO-A/MAO-B index of 27.7) was further subjected to reversibility and time-dependence inhibition studies, which disclosed a slow and irreversible inhibition of human MAO-B. Overall, the results support the suitability of the 4-triazolylalkyl propargylamine scaffold for exploring the design of multipotent anti-Alzheimer compounds endowed with irreversible MAO-B inhibitory activity.

New polycyclic dual inhibitors of the wild type and the V27A mutant M2 channel of the influenza A virus with unexpected binding mode.

Two new polycyclic scaffolds were synthesized and evaluated as anti-influenza A compounds. The 5-azapentacyclo[6.4.0.0(2,10).0(3,7).0(9,11)]dodecane derivatives were only active against the wild-type M2 channel in the low-micromolar range. However, some of the 14-azaheptacyclo[8.6.1.0(2,5).0(3,11).0(4,9).0(6,17).0(12,16)]heptadecane derivatives were dual inhibitors of the wild-type and the V27A mutant M2 channels. The antiviral activity of these molecules was confirmed by cell culture assays. Their binding mode was analysed through molecular dynamics simulations, which showed the existence of distinct binding modes in the wild type M2 channel and its V27A variant.

Assessing the suitability of the multilevel strategy for the conformational analysis of small ligands.

Predicting the conformational preferences of flexible compounds is a challenging problem in drug design, where the recognition between ligand and receptor is affected by the ability of the interacting partners to adopt a favorable conformation for the binding. To explore the conformational space of flexible ligands and to obtain the relative free energy of the conformation wells, we have recently reported a multilevel computational strategy that relies on the predominant-state approximation-where the conformational space is partitioned into distinct conformational wells-and combines a low-level method for sampling the conformational minima and high-level ab initio calculations for estimating their relative stability. In this study, we assess the performance of the multilevel strategy for predicting the conformational preferences of a series of structurally related phenylethylamines and streptomycin in aqueous solution. The charged nature of these compounds and the chemical complexity of streptomycin make them a challenging test for the multilevel approach. Furthermore, we explore the suitability of using a molecular mechanics approach as a source of approximate ensembles in the first stage of the multilevel strategy. The results support the reliability of the multilevel approach for obtaining an accurate conformational ensemble of small (bio)organic molecules in aqueous solution.

The complex of hypericin with ß-lactoglobulin has antimicrobial activity with potential applications in dairy industry.

Using a combination of molecular modeling and spectroscopic experiments, the naturally occurring, pharmacologically active hypericin compound is shown to form a stable complex with the dimeric form of ß-lactoglobulin (ß-LG). Binding is predicted to occur at the narrowest cleft found at the interface between monomers in the dimeric ß-LG. The complex is able to preserve the fluorescence and singlet oxygen photosensitizing properties of the dye. The equilibrium constant for hypericin binding has been determined as Ka=1.40±0.07µM(-1), equivalent to a dissociation constant, Kd=0.71±0.03µM. The complex is active against Staphylococcus aureus bacteria. Overall, the results are encouraging for pursuing the potential application of the complex between hypericin and ß-LG as a nanodevice with bactericidal properties for disinfection.

Tetrahydrobenzo[h][1,6]naphthyridine-6-chlorotacrine hybrids as a new family of anti-Alzheimer agents targeting ß-amyloid, tau, and cholinesterase pathologies.

Optimization of an essentially inactive 3,4-dihydro-2H-pyrano[3,2-c]quinoline carboxylic ester derivative as acetylcholinesterase (AChE) peripheral anionic site (PAS)-binding motif by double O → NH bioisosteric replacement, combined with molecular hybridization with the AChE catalytic anionic site (CAS) inhibitor 6-chlorotacrine and molecular dynamics-driven optimization of the length of the linker has resulted in the development of the trimethylene-linked 1,2,3,4-tetrahydrobenzo[h][1,6]naphthyridine-6-chlorotacrine hybrid 5a as a picomolar inhibitor of human AChE (hAChE). The tetra-, penta-, and octamethylene-linked homologues 5b-d have been also synthesized for comparison purposes, and found to retain the nanomolar hAChE inhibitory potency of the parent 6-chlorotacrine. Further biological profiling of hybrids 5a-d has shown that they are also potent inhibitors of human butyrylcholinesterase and moderately potent Aß42 and tau anti-aggregating agents, with IC50 values in the submicromolar and low micromolar range, respectively. Also, in vitro studies using an artificial membrane model have predicted a good brain permeability for hybrids 5a-d, and hence, their ability to reach their targets in the central nervous system. The multitarget profile of the novel hybrids makes them promising leads for developing anti-Alzheimer drug candidates with more balanced biological activities.

Easily accessible polycyclic amines that inhibit the wild-type and amantadine-resistant mutants of the M2 channel of influenza A virus.

Amantadine inhibits the M2 proton channel of influenza A virus, yet most of the currently circulating strains of the virus carry mutations in the M2 protein that render the virus amantadine-resistant. While most of the research on novel amantadine analogues has revolved around the synthesis of novel adamantane derivatives, we have recently found that other polycyclic scaffolds effectively block the M2 proton channel, including amantadine-resistant mutant channels. In this work, we have synthesized and characterized a series of pyrrolidine derivatives designed as analogues of amantadine. Inhibition of the wild-type M2 channel and the A/M2-S31N, A/M2-V27A, and A/M2-L26F mutant forms of the channel were measured in Xenopus oocytes using two-electrode voltage clamp assays. Most of the novel compounds inhibited the wild-type ion channel in the low micromolar range. Of note, two of the compounds inhibited the amantadine-resistant A/M2-V27A and A/M2-L26F mutant ion channels with submicromolar and low micromolar IC50, respectively. None of the compounds was found to inhibit the S31N mutant ion channel.